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Space-based experiments offer sensitivity to numerous unmeasured effects involving Lorentz and 
CPT violation. We provide a classification of clock sensitivities and present explicit expressions 
for time variations arising in such experiments from nonzero coefficients in the Lorentz- and CPT- 
violating Standard-Model Extension. 
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I. INTRODUCTION 

Unification of the fundamental forces in nature is ex- 
pected to occur at the Planck scale, rap ~ 10 19 GeV, 
where quantum physics and gravity meet. Performing 
experiments with energies at this scale is presently in- 
feasible, but suppressed signals might be detectable in 
exceptionally sensitive tests. Searching for violations of 
relativity that might occur at the Planck scale via the 
breaking of Lorentz and CPT symmetry is one promis- 
ing approach to uncovering Planck-scale physics [lj • 

At low energies relative to the Planck scale, observable 
effects of Lorentz violation are described by a general ef- 
fective quantum field theory constructed using the parti- 
cle fields in the Standard Model. This theory, called the 
Standard-Model Extension (SME) [2|, allows for general 
coordinate-independent violations of Lorentz symmetry. 
It provides a connection to the Planck scale through oper- 
ators of nonrenormalizable dimension [3| . CPT violation 
implies Lorentz violation T], so the SME also describes 
general effects from CPT violation. 

Various origins are possible for the Lorentz and CPT 
violation described by the SME. An elegant and generic 
mechanism is spontaneous Lorentz violation, originally 
proposed in the context of string theory and field theo- 
ries with gravity [j| and subsequently extended to include 
CPT violation in string theory 6] . Noncommutative field 
theories offer another popular field-theoretic context for 
Lorentz violation, in which realistic models form a subset 
of the SME Q • Lorentz violation has also been proposed 
as a feature of certain non-string approaches to quantum 
gravity, including loop quantum gravity and related mod- 
els of spacetime foam §j , the random dynamics approach 
H, and multiverse models [Toj . 

Various types of sensitive experiments can search for 
the low-energy signals predicted by the SME. In this 
work, we consider clock-comparison experiments with 
clocks co-located on a space platform, which are known 
to offer a broad range of options for Planck-sensitive tests 
of Lorentz and CPT symmetry . Promising pos- 

sibilities are offered by various experiments planned for 
flight on the International Space Statio n ( ISS), includ- 
ing the ACES [13, PARCS O, RACE and SUMO 
16] missions. The first three of these presently involve 



atomic clocks with 133 Cs, 87 Rb and a H maser ^tJ? while 
the fourth uses a superconducting microwave oscillator. 

Clock-comparison experiments in laboratories on the 
Earth [H [3 H3, HJ H3 have already demonstrated 
exceptional sensitivity to spacetime anisotropics at the 
Planck scale. These experiments monitor the frequency 
variations of a Zeeman hyperfine transition as the instan- 
taneous atomic inertial frame changes orientation. Typ- 
ically, a pair of clocks involving different atomic species 
and co-located in the laboratory is compared as the Earth 
rotates. Several other types of experiments are also sensi- 
tive to Planck-scale effects predicte d by the SME, includ- 
ing ones involving photon s IT2I l23l |24| |. hadrons pol |26| . 
muons 1^3, and electrons |28ll29||. 

In the present work, we perform a general analysis of 
clock-comparison experiments involving atomic clocks on 
a satellite such as the ISS. To take advantage of the rel- 
atively high velocities available in space, we incorporate 
leading-order relativistic effects arising from clock boosts. 
A framework for general calculations of this type is pre- 
sented, and detailed expressions that allow for satellite 
and Earth boosts are derived for observables in a stan- 
dard satellite mode. Estimates are provided of the sensi- 
tivities of experiments attainable on the ISS. 

The paper is organized as follows. Section [H] consid- 
ers some aspects of the frequency shifts due to Lorentz 
violation that are experienced by a clock in a single in- 
ertial frame. In Sec. 11111 we establish the link between a 
noninertial clock frame on a space platform and the stan- 
dard Sun-based frame. Section llVI presents methods for 
extracting measurements of coefficients for Lorentz vio- 
lation from experimental data and estimates sensitivities 
for ISS-type missions. We summarize in Sec.[V] Details 
of some calculations are provided in some appendices. 
Throughout this work, we adopt the notation of Refs. 

u na- 



il. BASICS 

Any Zeeman transition frequency u> used to study 
Lorentz and CPT violation can be written in the form 
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Here, B 3 is the magnitude of the external magnetic field 
when projected along the quantization axis, /(S 3 ) is the 
transition frequency according to conventional physics, 
and Slo contains all contributions from Lorentz and CPT 
violation. All orientation dependence is contained in B 3 
and Slo; in particular, the function / has no orientation 
dependence except through B3. Typically, / depends on 
magnetic moments, angular-momentum quantum num- 
bers, and similar quantities. For definiteness in what 
follows, we suppose / is invertible in a neighborhood of 
the magnetic fields of interest [3(|, and denote the in- 
verse of / by / _1 (x). Also, we work at all orders in B3 
but neglect effects of size o(B 3 5lo) and o(Slo 2 ), which are 
known to be small. 

For the transition (F,m F ) — > (F',m' F ), the frequency 
shift Slo can be written as 

Slo = SE{F,m F ) - 6E(F',m' F ), (2) 

where the atomic energy shifts SE(F, m F ) are induced by 
Lorentz and CPT violation. These shifts can be calcu- 
lated directly within the SME using standard perturba- 
tion theory, by obtaining the individual energy shifts for 
each constituent particle and combining the results. In 
the clock frame, they are determined at leading order by 
a few combinations of SME coefficients for Lorentz viola- 
tion, conventionally denoted as bf, df, g™, c%, g™, where 
the superscript w is p for the proton, n for the neutron, 
and e for the electron. These are the only quantities in 
the clock frame that can in principle be probed in clock- 
comparison experiments with ordinary matter |22|. 

In the clock frame, the atomic energy shift for state 
\F,m F ) can be written 

6E(F,m F ) = m F ^(W + W^Mj) 

W 

+ m F Y,{l w ~c w q +\ w ~g™) ■ (3) 

w 

Here, fh F and m F are specific ratios of Clebsch-Gordan 
coefficients, while /3 W , 5 W , k w , j w , X w are specific expec- 
tation values of combinations of spin and momentum op- 
erators in the extremal states \F,m F = F). For present 
purposes, the details of these quantities are unnecessary; 
they are given in Eqs. (7), (9), (10) of Ref. H^. 

Clock-comparison experiments typically involve two 
clocks and corresponding transitions A, B with frequency 



shifts Sloa, Slob, located in an external magnetic field B 3 . 
The experimental signal of interest is a modified differ- 
ence between frequency shifts of the form Sloa — vSlob, 
where v is an experiment-specific constant related to the 
gyromagnetic ratios of the two clocks. In typical arrange- 
ments, v is such that this signal vanishes in the absence 
of Lorentz violation. Note that the two transitions may 
involve the same atomic species. 

To bridge experiment and theory, it is useful to intro- 
duce a modified frequency difference lo^ that represents 
the signal for a large class of experimental situations and 
offers a direct link to coefficients for Lorentz violation in 
the SME. For the two clock transitions A, B with fre- 
quencies to a, u>b written in the form (JJJ, define to* by 

J :=u A ^ f A {f B 1 {u JB )). (4) 

By construction, ur vanishes in the absence of Lorentz 
violation. To the order in Slo at which we work, this 
implies lo^ is independent of the external magnetic field 
even in the presence of Lorentz violation. It is therefore 
reasonable to adopt this definition of to* as the ideal ob- 
servable for Lorentz and CPT violation. In what follows, 
we first obtain a general theoretical expression for ufi and 
then consider some experimental issues. 

We next show that w" is determined theoretically by 
the equation 

u>" = Slo a — vSlob, (5) 

where 
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and Slo a, Slo b are given by Eq. J2J. This expression for v 
is valid to all orders in B3. When combined with Eqs. Q), 
@, and ©, the above two equations allow calculation 
of aA 

To prove Eqs. (J5J and ©, we proceed as follows. For 
each transition of the form (|TJ, define an effective mag- 
netic field B cS = / _1 (o;). In the special case of no 
Lorentz violation, B cfi is identical to the actual mag- 
netic field S3, so the difference Bf - Bf between the 
transitions A, B is zero. However, in general we have 



Bf - Bf = f^(Lo A ) - f B \oo B ) = Fa 1 [/a (-S3) + Sloa] - fs 1 [/s (S3) + 5lo b ] 



= Sua 



dx 



— Slo 1 



dfs 1 



x=fA(B 3 ) 



dx 



o(Slo) 2 



(7) 



x=f B {B 3 ) 



where Taylor expansions in Sloa and Slob have been performed. This implies 
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where another Taylor expansion has been performed. 
Within factors of size o(B 3 5lo), we can set B3 = on 
the right-hand side of this equation, except for the term 
/aI/b^^-b)]- Applying the identity (df' 1 /dx)\ x=f{B3 ) = 
(df/dB 3 )\ B l then yields Eqs. © and ©. 

As a first example of calculation with these results, 
consider the special case of linear dependence on B 3 . 
Suppose for each transition we can write f(B 3 ) — c+pB 3 , 
where c and p are constants for each transition. Then, 
we find 
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In this case, it suffices to study the combination 
lua — Pa^b/pb and neglect the constants, since clock- 
comparison experiments are only sensitive to orientation- 
dependent effects. 

For a more complicated example, consider the special 
case of a quadratic dependence on B 3 . Suppose for each 
transition we can write f(B) = c+ pB + pB 2 , where again 
c, p, p are constants for each transition. As always, lo^ is 
relatively simple when expressed in terms of frequency 
shifts for Lorentz violation: u* — Sloa — PaSlob/ Pb- 
However, in terms of the individual frequencies or is 
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constant terms. 



P 2 b 



4p B (^B - c B ) 



(10) 



Note that the previous linear example is a nontrivial limit 
of this one because / _1 behaves badly as p — ► 0. 

A clock-comparison experiment to probe Lorentz vi- 
olation can proceed in several ways. The most direct 
method is to measure loa and ub at each instant. The 
results are then combined according to Eq. 10} to give an 
experimental value of w", which may be compared to the 
theoretical calculation in Eq. (J5J. A potentially signif- 
icant disadvantage of this method is that achieving the 
desired sensitivity requires exquisitely precise knowledge 
of the functions /a and fg and the parameters on which 
they depend. 

A different procedure can be adopted that requires no 
knowledge of the functions Ja and fs- Suppose u>b is 
forced to be constant, perhaps by applying a feedback 
magnetic field [l!4|2(|. Then, /^[/^(wb)] is constant, 
so ct>" = loa up to a constant irrelevant for experimental 
purposes. Thus, if lob is held constant and the transitions 
A and B involve clocks subject to the same instantaneous 
magnetic field, it follows that loa — lo^ = Sloa — vSlob- 
Then, loa is sensitive purely to Lorentz-violating effects 




FIG. 1: Orbit of Earth in Sun-based frame 



and can be interpreted without detailed knowledge of Ja 
and fs- This procedure may offer practical advantages 
for experiments in environments with fluctuating mag- 
netic fields such as those anticipated for the ISS experi- 
ments. 



III. FRAME TRANSFORMATIONS 

In a clock-comparison experiment, the instantaneous 
clock frame is continuously changing due to the orbital 
and rotational motion of the space-based laboratory [3l| . 
The quantities bf, df, g™, , g™, therefore vary in 
time, with frequencies determined by the orbital and 
rotation periods of the laboratory. This time varia- 
tion can be obtained explicitly by converting 63, df , 
9d j Cq-> 9qi from the laboratory frame with coordinates 
(0, 1,2,3) to a specified nonrotating frame with coordi- 
nates (T,X,Y,Z). 

Following Refs. 0,0], in this work we adopt for the 
standard nonrotating frame a natural Sun-centered ce- 
lestial equatorial frame. This frame is approximately in- 
ertial over thousands of years. It is therefore suitable 
for the study of leading-order boost effects due to the 
Earth and satellite orbital motions. The results of all 
clock-comparison experiments to date can be regarded 
as having been reported in this frame. 

In the Sun-based frame, the spatial origin coincides 
with the center of the Sun. The unit vector Z is paral- 
lel to the Earth's rotational axis, X points to the ver- 
nal equinox on the celestial sphere, and Y completes the 
right-handed system. The time T is measured by a clock 
fixed at the origin, with T = chosen as the vernal 
equinox in the year 2000. Note that the vectors X, Y 
lie in the Earth's equatorial plane, which itself is at an 
angle of 77 « 23° to the Earth's orbital plane. Note also 
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FIG. 2: Parameters for definition of satellite orbit. To sim- 
plify the presentation, Earth is pictured as if it were translated 
to the Sun-frame coordinate origin. 



that the Earth is on the negative X axis at time T = 0. 
See Fig. [2 

For a space-based experiment, the time variation of the 
clock frequency is determined by the satellite orbital and 
rotational motions. To extract the leading-order effects 
relevant for experiments on the Earth and on the ISS, it 
suffices to approximate the orbits as circles. Any elliptic- 
ity introduces time dependence at higher harmonics of or- 
bital frequencies, suppressed by even powers of the orbit 
eccentricity e 2 . For example, a time dependence propor- 
tional to coswi under the circular-orbit approximation 
generates an order-e 2 dependence ~ e 2 cos 3u>t for an el- 
liptical orbit. These harmonics appear only at subleading 
order for any quantity that they modify. For present pur- 
poses, the circular approximation is reasonable because 
e% ~ 0.029 for the Earth's orbit and e 2 ~ 0.032 for the 
ISS orbit. However, dedicated satellite missions could 
have strongly elliptical orbits, in which case the higher 
harmonics would be of interest. 

Under the circular approximation, the parameters of 
the Earth's orbit are the mean orbital radius ii® and the 
mean orbital angular frequency CIq. The mean Earth or- 
bital speed is /3® = R&Q^. The parameters for a circular 
satellite orbit around the Earth are taken as the mean 
orbital radius r s , the mean satellite orbital angular fre- 
quency uj s , the angle £ between the Earth's rotation axis 
Z and the satellite orbital axis, the azimuthal angle a 
between the satellite and the Earth orbital planes, and a 
conveniently chosen reference time To at which the satel- 
lite crosses the equatorial plane on an ascending orbit. 
See Figure |21 It is also useful to introduce the satellite 
time measured in the Sun-based frame, T s = T — To. 
Note that the mean satellite speed with respect to the 
Earth's center is (3 S — r s oj s . For special limiting orbits, 
uj s reduces to the usual sidereal frequency [33 . Note also 
that various perturbations typically cause a to precess. 

The rotational motion of the satellite is specified by 
giving its orientation as a function of time. Two flight 
modes are commonly considered [33j , often denoted XVV 



and XPOP. In XPOP mode, the satellite orientation is 
fixed in the Sun-based frame as it orbits the Earth. All 
clock signals from Lorentz violation are due to boosts as- 
sociated with the satellite orbital motion in this frame, 
so they are suppressed by at least one power of /3® or (3 S . 
In contrast, for the XVV ("airplane") mode, the satel- 
lite rotates once in the Sun frame each time it orbits the 
Earth, so its orientation is fixed relative to the instanta- 
neous tangent to the satellite's circular orbit about the 
Earth. Clock signals in this mode are due to both rota- 
tions and boosts, so they are sensitive to a wide variety 
of Lorentz-violating effects. In what follows, we focus on 
the XVV mode. 

In the space-based laboratory, the coordinate system is 
defined as follows 0, 0] . The 3 axis is taken along the 
satellite velocity with respect to the Earth. The I axis is 
chosen to point towards the center of Earth. The 2 axis 
completes the right-handed system and is oriented along 
the satellite orbital angular momentum with respect to 
the Earth. The clock orientation in the laboratory is 
typically determined by an applied magnetic field, which 
establishes a quantization axis. For definiteness, we take 
the quantization axis as the 3 axis in this work. Other 
choices of quantization axis can readily be calculated by 
our methods [l^. Although the detailed time- varying 
signals are different, no additional sensitivities to Lorentz 
violation are obtained with other choices. 

Combining information from the above frame, mode, 
and orientation choices permits the construction of the 
explicit transformation T between the Sun-based and 
laboratory frames. Acting on vector components, the 
transformation can be regarded as a matrix with compo- 
nents T^" that depend on various velocities, frequencies, 
angles, and Sun-frame times. The derivation of this ma- 
trix is provided in Appendix ^ With this matrix, the 
explicit time dependence of the quantities 63 , , g"%, c™ , 
g™, can readily be calculated in terms of the Sun- frame 
coefficients <zg, 6g, <f n , dg n , e%,f%, gf ns , H£ n appear- 
ing in the fermion sector of the SME, where 2, II, X are 
indices spanning the Sun-frame coordinates (T, X, Y, Z). 
For example, 63 = T^bs, and do3 — Tq^T^ 11 dsn- 
Due to the relatively involved spiral nature of the satel- 
lite trajectory as observed in the Sun-frame, the resulting 
explicit expressions for the quantities bf, dg , g™, c™, g™, 
are somewhat lengthy. It turns out to be simpler and nat- 
ural to express these in terms of certain special "tilde" 
combinations of Sun-frame coefficients for Lorentz vio- 
lation (34|. These combinations are listed in Appendix 
IbI For each of the three species, 40 independent Sun- 
frame tilde coefficients play a role at the level of zeroth- 
and first-order relativistic effects considered here. There 
are therefore 120 linearly independent degrees of freedom 
that can be probed in clock-comparison experiments with 
ordinary matter at this relativistic order. Note that for 
each species the SME coefficients a^, b^, s d^, e M , / M , 
gx^v , H)lv contain a total of 44 physically observable co- 
efficients at leading order in Lorentz violation once un- 
physical field redefinitions have been fixed 0, |3f|, |3(J , so 
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TABLE I: Dependence of clock-frame coefficients on satellite time T a and on Sun-frame tilde coefficients. 



four additional Sun-frame tilde coefficients are required 
to form a complete set of physical observables for clock- 
comparison experiments. However, these can appear at 
most as subleading-order relativistic effects with signals 
suppressed by two powers of the velocities /3®, (3 S and are 
therefore not considered in this work. 

The resulting expressions for b™ , dg , g™, c™ , g™, are 
given in Appendix |0 Each equation is a linear combi- 
nation of Sun-frame tilde coefficients. The multiplicative 
factors are constants of order f , sines or cosines of the 
angles a, 2a, £, 2£, 77, and time oscillations involving 
sines or cosines of w s T s , 2ui s T s , f2®T. Note that the 
terms involving f2® vary relatively slowly with time be- 
cause u s 3> Qq. The same is true of any precession time 



dependence in the orbital angle a. Note also that the 
usual nonrelativistic dependence \22\ is recovered in the 
nonrclativistic limit /3® — > 0,/3 s — * 0. 

Insight into the content of these equations can be ob- 
tained by separating each according to distinct satellite- 
frequency dependences and classifying the resulting 
terms according to velocity dependence. Since (3® ~ 
I(T 4 for the Earth and (3 S ~ I(T 5 for the ISS, the 
terms linear in the velocities are suppressed relative to 
the zeroth-order ones. Table U lists the decomposition 
of the equations in Appendix [C] in accordance with this 
scheme. As an explicit example, consider the variation 
of c q with the fundamental satellite frequency oj s . This 
is contained in the full expression for c q presented in Eq. 
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(|C4|) , from which the relevant terms can be extracted and 
rearranged in the form 

c q D P s (2s a c c c TX - 2c a c c c T Y ~ 2s c c tz )cosuj s T s 

+(3 s (2c a c T x + 2s Q c T y) sinuj s T s . (11) 

Table [I] separates the sine and cosine dependence of 
this expression and lists factors of (3 S in the appropriate 
columns for the coefficients cpx, cpy, ctz- In general, 
this type of structural information about the equations 
in Appendix El is useful in establishing sensitivities for 
different experiments. 



IV. SIGNALS AND SENSITIVITIES 

At this stage, we can use the time dependence of the 
quantities 63 , d™ , g% , c™, g™, derived in the previous sec- 
tion to study the signals in clock-comparison experiments 
involving various atomic transitions. We focus specifi- 
cally on transitions (F,mp) — > (F 1 , m' F ) in species sched- 
uled for flight on the ISS: 87 Rb, 133 Cs, and H. These have 
an even number of neutrons and total electronic angular 
momentum J = 1/2. The generalization of our results to 
nuclei with an odd number of neutrons is straightforward. 

The Lorentz-violating contribution Sto to the frequency 
of the transition (F,mp) — * (F',m F ) is given by Eqs. 
|j2J and ©■ With the above assumptions, this frequency 
shift can be expressed as: 

Sw = s? P p (lj)bl + 5 p (lj)d% + n p {lj)g p d 

+4 [/3 e (0i/ 2 )^ + <5 e (0 1/2 )dl + Ke (0 1/2 )d .(12) 



The second case has F = F' 

P t~- i>l \ (7+l)(2I-l 

s l = ( m F - m F y J(2 ; ) +1) 
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(J+l)(2/-l 



|, and we find 
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1(21+1) 



1(21+1) 

Amp 



1(21+1) 

m'p 



' F' 



7(2/+l)(2/-l) 



mp 



1-2/ 



l+2I 77l F 

Amp + 



21+1 

The final case has F = I 

(J+l)(27-l) 



27+1 

2 ' r 



(15) 



1(21+1) 



mp 



Amp 



21+1 

(J-l)(2J+3) 
/(2/+1) 

3 



mp 



1(21+1) 

m'p 



i, for which 



m F , 



1(21+1 

1-2/ 



A? 



12 



/(2/+l)(2/-l) 



1 p' 



1+21 
2 



mp 



m 1 



21+1 



Amp 



21+1 



mp. 



(16) 



Various results can be obtained from these expressions 
and Eq. (|12|l . For example, it follows directly that a 
nonzero signal occurs for all AF = ±1, Am^ = tran- 



In this equation, each lj refers to the Schmidt nucleon. sitions except for the special case m F 



Except for the quantities sj outside the brackets, all vari 
ables in Eq. (|12fl are those appearing in Eq. The 
specific values of the quantities j3 w , q^, S w , k w , X w are 
given as Eqs. (11) and (12) of Ref. [22]. The values of 
depend on the transition, and formulae for them are 
given below. Note that similar equations valid for more 
general atoms would also involve s|, s™ , and terms. 

The expressions for the sj can be classified according 
to the possible values of F and F' . There are four cases 
of interest. For each case, we give the expressions first in 
terms of combinations of Clebsch-Gordan coefficients and 
angular-momentum quantum numbers, and then directly 
in terms of m F and m' F . In all cases, we define Amp :— 
m F — m' F and Amp := mp — (m' F ) 2 . 



i' F = 0. This 

demonstrates that the standard clock transitions are in- 
sensitive to Lorentz-violating effects, in agreement with 
previous results j2^. Useful special cases of immedi- 
ate relevance to experiments on the ISS can also be ex- 



tracted. Thus, for a 133 Cs clock with I 
we find: 



F 



3, 



s 2 



= jAmp 

= 2% Am F 

= \Amp 



jj riip, 
- Jiim'p) 2 



Similarly, for a 87 Rb clock with / 
obtain: 



The first case has F = F' = I - 



p — 



p — 



m p 
friF 
rhp 



op 

1/ 



27+1 
3 



for which we obtain 



Amp, 



= ^Amp 
= \ Am| 
= \Amp 



1 F ■ 



F 



(17) 
1, we 



(m' F ) 
m'p. 



(18) 



7(27+1) 



Am 2 



F- 



IS 



27+1 



Amp. 



(13) 



For an H clock or maser, s\ — s\ = 1. However 
irrelevant: the proton is in an I = 1/2 state, so there is 
no quadrupole effect and the quantities j p and A p vanish. 
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133 Cs 


133 Cs 


133 Cs 


133 Cs 


133 Cs 


87 Rb 


87 Rb 


87 Rb 


87 Rb 


87 Rb 


X H 


transition 


(4,0) 


(4,1) 


(4,-1) 


(4,0) 


(4,0) 


(2,1) 


(2,-1) 


(2,1) 


(2,0) 


(2,0) 


(1,1) 




-(3,0) 


- (3, 1) 


-(3,-1) 


-(3,1) 


-(3,-1) 


-(1,1) 


-(1,-1) 


-(1,0) 


-(1,1) 


-(1,-1) 


-(1,0) 


I 


7/2 


7/2 


7/2 


7/2 


7/2 


3/2 


3/2 


3/2 


3/2 


3/2 


1/2 


Z 


55 


55 


55 


55 


55 


37 


37 


37 


37 


37 


1 


N 


78 


78 


78 


78 


78 


50 


50 


50 


50 


50 





Schmidt 


97/2 


97/2 


97/2 


97/2 


97/2 


PZ/2 


V-i/2 


PZ/2 


PS/2 


P3/2 


Sl/2 


nucleon 
























e~ state 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Sl/2 


Pp 


[§] 






[f] 


[|] 


hi] 


hi] 


hi] 


hi] 


hi] 


-1 


Iv 


[-£tfp] 




hi*,] 


H K p\ 






1-TsKp] 




[-js k p] 


1-TsKp] 





s p 

fvp 




[~hKp] 


[f 9 K P ] 




[§*>] 


[W 
HKp] 


[\K P ] 
H K p\ 


[\K P ] 


\\Kp\ 
H K p] 


HKp] 


S*p 
-\K P 


Xp 


[0] 


[0] 


[0] 


[0] 


[0] 


[0] 


[0] 


[0] 


[0] 


[0] 





Pe 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


7e 



































8e 


\K e 


\K e 


\K e 


\K e 


\K e 


\K e 


\K e 


\K e 




\K e 


\K e 


fi>e 




-\K e 


-\K e 


-\K e 


-\K e 


-\K a 




-\K. 


-\K e 


-\K e 


-\K e 


Xe 



































-? 





-1/14 


1/14 


-9/28 


9/28 


-1/3 


1/3 


1/2 


-5/6 


5/6 


1 







-1/14 


-1/14 


-5/28 


-5/28 


-1 


-1 


1/2 


-3/2 


-3/2 




Si 





1/2 


-1/2 


1/4 


-1/4 


1 


-1 


1/2 


1/2 


-1/2 


1 



TABLE II: Parameters for transition frequencies for experiments with Cs, Rb, and H clocks. 



Table [H] summarizes some useful results for species 
scheduled for flight on the ISS. The first few rows of this 
table identify the transition and list various properties of 
the species involved. The nuclear spin is /, the proton 
number is Z ', and the neutron number is N . The follow- 
ing entry fixes the proton determining the ground-state 
properties of the nucleus following the nuclear Schmidt 
model |3tJ, together with its associated orbital and to- 
tal angular momentum. The electronic configuration is 
also given. Ten rows list the relevant parameters /3 W , 
8 W , K w, 7t«, Atu, with values in brackets obtained un- 
der the assumptions of the Schmidt model. We define 
K p = {p 2 )/m,p, which can be regarded as twice the ki- 
netic energy per mass of the Schmidt-model proton, and 
define K e similarly for the valence electron. An estimate 
gives K p ~ 10~ 2 for all species except 1 H, for which 
K p ~ 10~ n , and K e ~ 10~ 5 . Finally, the numerical 



values of the sj are listed. Where these are nonzero, 
a clock-comparison experiment with the specified transi- 
tion is sensitive to Lorentz violation. 

At this stage, enough information is at hand to extract 
estimated experimental sensitivities. Suppose the results 
of an experiment measuring the modified frequency dif- 
ference ur of Eq. (0} are fitted to the form 

or = constant -I- 27re li x coslo s T s + 27re 1 y sinu) s T s 

+2ne 2 ,x cos2uj s T s + 2ks 2 .y sin 2uj s T s . (19) 

Nonzero values of any of the e a ,j indicate Lorentz vi- 
olation. Denote by e a the minimum of {|£ ,x|, |£a,r|}- 
Then, combining the theoretical analysis above yields the 
following predicted dependence of e a on coefficients for 
Lorentz violation, atomic and nuclear parameters, and 
geometrical factors: 



27T£i = 



E 



W A fit - VS? B pZ Wn + ^ [ s ^ {fj A + S A + K A } vsr B {p B + § B + ^ (& ™ 



[sr A (tt + s A ) - + O] d%, ay K , Hy T ) + K A s A - v S y B s*)(dy) 

+ (^ A ^t - vs? B K*)(g%j) + [P m sT A (S A + ni) - /3 ffi ^r B (5f + «£) + Ps(^ A \i ~ vs^^)] {g w j K ) 



+Pe L s i \Pw + K w) - vs l [Pw + K w)\ (9c ) + Ps(S 2 ~/ w - VS 2 ~f w )(C TJ ) 



(20) 
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2tT£ 2 = 



E 



[p s s? A (/3 A + 6 A + 4) - Psvs^i^ + S B + k b ) + /3 @ (s% A \i - vs™ B \ B )} (b%,~9%) 



+f3 s [s? A (j3 A + S A ) - v S ? B (^ + 5 B )] {d w _,d w Q ,d w JK ) + (3 s (si A p A - vsr B p B )(Hy T ) 
+ [P s s? A (6 A + k A ) - P s vsy B (5 B + k b ) + Ms^Xi - vs™ B \ B )] (gy K ) 



+ [P S S? A (P A + 4) - (3sVsr B (P B + n B )+f3 @ ( S ^ A K 



vs% B ' B 



wA„A 



wB B\/ ~w -w ~w\ i / w/i. \ j\ w ts \ &\ t ~w ~w ~w \ 

1w)k c -> c qi c j) + ( s 2 K~vs 2 X w ){g-,gQ,9Tj) 



wA \A 



wB\B\{~w ~w ~w 



(21) 



In these equations, superscripts A and B indicate quanti- 
ties evaluated for transitions A and B, respectively, and 
J takes the values (X, Y, Z). The coefficients for Lorentz 
violation enter as somewhat lengthy linear combinations 
of the type appearing in the equations of Appendix [U] 
These explicit combinations are omitted here for brevity, 
being replaced instead with parentheses containing only 
the specific coefficients for Lorentz violation involved. 

The above form of the equations is useful despite the 
brevity because it allows relatively straightforward con- 
sideration of sensitivities to the Sun-frame tilde coeffi- 
cients for Lorentz violation. We adopt here the strat- 
egy of Ref. [23, m which numerical sensitivities are ob- 
tained within the Schmidt model under the plausible as- 
sumption that no substantial cancellations occur among 
contributions from different Sun-frame tilde coefficients 
for Lorentz violation. For example, if e\ is an experi- 
mental sensitivity to the time variation of u>\ then Eq. 
(|20|l implies the experiment has sensitivity to each \b u j \ of 
- 2tt£ 1 (s5< m /3^ - vsf 3 ^)- 1 . Similarly, the sensitivity 
to \V£\ is ~ 2ne 1 ^ 1 {s? A ((3 A + S A + k a ) - vsf B (0 B + 
S B + and so on. To obtain crude order of mag- 

nitude numerical estimates, it suffices to approximate 
f3 s ~ 10~ 5 , /3® ~ 10~ 4 , and to estimate nonzero values of 
the other parameters as follows: (3 W ~ 1, s™ ~ 1 for all 
species; S, k, 7, A ~ 10~ 2 for protons except for 1 H where 
the nonzero values are only 6, k ~ 10~ n for the proton; 
and 5, A ~ 10~ 5 for electrons. Sensitivity estimates 
of this type are reasonable provided the various angles a, 
2a, C, 2£, 77, ri©T lie away from multiples of ir/4. The 
orientation of the quantization axis within the satellite 
then makes little difference to the sensitivity. However, 
for any angles near a multiple of 7f/4, sensitivity to one 
or more of the Sun- frame tilde coefficients can be lost. 

Table IIIII lists estimated sensitivities to Sun-frame 
tilde coefficients for Lorentz violation that might be at- 
tained in the planned space-based clock-comparison ex- 
periments with 133 Cs and 87 Rb clocks. The base-10 log- 
arithm of the sensitivity per GeV is shown for each coef- 
ficient for Lorentz violation and for each particle species. 
For definiteness, the clock sensitivity has been taken as 
£1,2 ~ 50 /xHz, which is comparable to that attained in 
a ground-based experiment with 133 Cs 0, but the re- 
sults shown are readily scaled for other values of £1,2- 
A star in the table indicates a combination for which 
there is no sensitivity according to the nuclear Schmidt 



Coefficient Proton Neutron Electron 



bx, by 


-27[-27] 


[-31] 


-27[-29] 


bz 


-27 


- 


-27[-28] 


br 


-23 


- 


-23 


in 


-23 


- 


-23 


Hjt 


-23 


- 


-23 


d± 


-23 


- 


-23 


dQ 


-23 


- 


-23 


djK 


-23 


- 


-23 


d~x, d Y 


-25 [-25] 


[-29] 


-22[-22] 


dz 


-25 




-22 


9DX,gDY 


-25 [-25] 


[-29] 


-22[-22] 


9DZ 


-25 




-22 


9JK 


-21 




-18 


9 c 


-23 




-23 


ctj 


-20 






C- 


-25 


[-27] 




CQ 


-25 






cx, c Y 


-25 


[-25] 




cz 


-25 


[-27] 




ctj 


-21 






g- 


*W 


M 




§q 


* 






grx, 9ty 




W 




grz 


*W 


[*] 





TABLE III: Estimated sensitivity to coefficients for Lorentz 
violation for ISS ex per iments with 133 Cs and 87 Rb clocks. 
Existing bounds (lSi Hfll I2H l2ll l2flj are shown in brackets. 



model but probable sensitivity in a more realistic nuclear 
model. A value in brackets indicates an existing bound 
from an Earth-based experiment [Tl HI I M I2I |29 | . 
Given the approximations described above, some caution 
in interpretation of the details of this table is advisable. 
Nonetheless, the table provides a measure of the broad 
scope of space-based tests of Lorentz symmetry, and it 
shows that Planck-scale sensitivity for a wide spectrum 
of relativity tests is attainable. 

Space-based experiments exploring the photon sector 
of the SME have been studied elsewhere 12], including 
the SUMO experiment with superconducting microwave- 
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cavity oscillators that is presently scheduled for flight. 
Experiments of this type could be profitably combined 
with the clock-comparison experiments discussed in the 
present work. For example, at the time of writing PARCS 
and SUMO are planned for simultaneous flight. Several 
configurations of interest could then be considered, in- 
cluding operating PARCS on a Lorentz-insensitive line 
as a reference for relativity tests with SUMO, or seek- 
ing uj s and 2oj s signals in a configuration with both the 
atomic clock and the cavity operating in modes sensitive 
to Lorcntz violation. 



V. SUMMARY 

This work has studied clock-comparison experiments 
on a space-based platform, with specific emphasis 
placed on forthcoming experiments on the International 
Space Station. The theoretical framework adopted is 
the Standard-Model Extension, which describes general 
Lorentz and CPT violation. The analysis yields predic- 
tions for signals at the ISS orbital and double-orbital fre- 
quencies, along with slower variations associated with the 
Earth orbital motion. 

The formalism we have presented applies to any space- 
based experiment with atomic clocks and incorporates 
relativistic effects at first boost order. We have derived 
explicit expressions for the observable effects in the spe- 
cial cases of 133 Cs, 87 Rb, and H clocks on the ISS, which 
are currently planned for flight in the PARCS, ACES, and 
RACE missions. These results, which involve the fermion 
sector of the SME, complement the photon-sector anal- 
ysis of Lorentz-violation sensitivity performed for the 
planned SUMO experiment with microwaves on the ISS. 

We have obtained estimates for the attainable sensi- 
tivities with these atomic-clock missions, listed in Ta- 
blc IIIII Numerous currently unmeasured coefficients for 
Lorentz violation could be studied in these experiments. 
The results demonstrate that experiments of this type 
offer potential sensitivity to violations of relativity with 
Planck-scale reach. 
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APPENDIX A: TRANSFORMATION FROM 
SUN-BASED TO SATELLITE FRAME 

In this appendix, we derive the transformation matrix 
introduced in Sec. IIIII that maps Sun- frame quanti- 
ties to laboratory- frame ones. The transformation T can 
be expressed as the composition of a boost A from the 
Sun frame to the (nonrotating) rest frame of the center 
of the satellite, followed by a rotation 1Z from the (non- 
rotating) rest frame of the center of the satellite to the 
(rotating) lab frame. Each constituent transformation 
depends on time, and each is understood to be instan- 
taneous. We first obtain an expression for the satellite 
position in the Sun frame, then use this to derive the 
instantaneous satellite velocity in the Sun frame, and fi- 
nally combine information to construct the desired over- 
all transformation. The conventions we adopt are those 
given in Refs. [HQ3. 

The position X@ of the center of the Earth in the Sun- 
based frame is given by 



X, 



f -i?©cosf2©T 

— i?^ cos X] sin 0©T 
y — i?© sinTysinfi^T 



(Al) 



The satellite position X s in this frame is obtained by 
adding the position of the satellite with respect to the 
Earth, which gives 



— cos ri©T + r s cosacosid s T s — r s cos ( sin a sin uj s T s 
X x = | — i?© cos rj sin f2©T + r s sin a cos uj s T s + r s cos a cos £ sin u s T s 
— i?0 sinTysinf^T + r s sin£sinw s T s 



(A2) 



Disregarding rotations for the moment, the boost A 
from the Sun-based frame to the nonrotating instanta- 
neous satellite rest frame is determined by the velocity 



V = dX s /dT. To lowest order in \V\, this is 



/ 1 



A 



V x 
V z 

V v z 



Vx Vy Vz \ 





1 
1 / 



(A3) 
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The required rotation 1Z from the nonrotating instan- 
taneous satellite rest frame to the laboratory frame on 
the satellite may be calculated using the velocity and 
acceleration of the satellite with respect to the Earth, 
d(X s -X®)/dT ~ z and d 2 (X s ~ X^/dT 2 ~ x, and the 
requirement that 



T Y 
■to 

7i T 



= -f3 e c n cnT + (3 s c a C(; cos lo s T s - (3 s s a sinw s T s , 
= -P m s v cor + PsS,; COSUJ s T s , 

= (SaCjjCftT - C a SQ, T )l3® COSUJ s T s 











— —c a cosuj s T s 


f s Q c c sinw s T s , 










— —s a cosu) s T s 


- c a q sinw s T s , 




-« 


[?)■ 


(A4) r x z 


= -s c sinw s T s , 










J-2 







Given A and 1Z, the overall instantaneous transforma- 
tion T from the Sun-based frame to the laboratory frame 
can be found: 



r 2 



X 



T = 



(1 





Vo 



o o \ 



r 2 z 



K 



A 



(A5) 



With the approximations of Sec. lIIII the components 
of the transformation matrix (|A5() are found to be 

% T = 1, 



= c o 

= /3 S + (s a c n cn T - c a so T )/3e sinw s T s 
+ (-c q c i( c c cot - SQSnC^T - s a C£SnT)0@ cosuj s T s , 
= —c a smuj s T s — s a c^ cos u) s T s , 
= c Q c c cosuj s T s - s a sinw s T s , 
= s^coslu s T s . (A6) 



x 



/3esnr - /3 s s a c c cos uj s T s - f3 s c a sinw s T s 



In these equations, the abbreviations s x = sin a;, c x 
cos x and QT = £1®T are used. 



APPENDIX B: SUN-FRAME COEFFICIENTS 

The Sun-frame tilde coefficients are defined as follows: 

bj = bj - \ejklHkl - m(d JT - \ejkl9klt), b T = b T + mg X YZ, 

g T = b T - m(g X YZ - 9yzx - 9zxy), 

Hxt = Hxt + iri(dzY — 9xtt — 9xyy), Hyt — Hyt + m(dxz — 9ytt — 9yzz)i 

Hzt = Hzt + m(dyx — gzTT — gzxx), 

d± = m(d X x ±rfyy), d Q = m(d X x + d Y v - 2d Z z - 9yzx - 9zxy + 2g X vz), 

dj = m(d TJ + \d JT ) - jEjklHkl, d~YZ = m(d Y z + d Z Y - 9xyy + 9xzz), 

dzx = m(d Z x + d X z ~ 9yzz + 9yxx), dxY = m(dxY + d Y x — 9zxx + 9zyy), 

9c = m(g X YZ - 9zxy), 9- = m(g X rx ~ 9yty), 9q = m(g X TX + 9yty - 2gzrz), 
9tj = m \ejkl\9ktl, 9dj = -bj + m£ JKL {g KTL + ^klt), 
9.JK = rn(g JT T + 9jkk) (no K sum, J ^ K), c Q = m(c x x + c YY - 2c Z z), 
c_ = m(c X x - c Y y), cj = m\ejKL\cKL, c T j = m(c T j + c JT ). (Bl) 

Indices J, K, L run over Sun-frame spatial coordinates X, Y, Z. The usual summation convention holds except where 
indicated. The totally antisymmetric tensor sjkl is defined with Exyz — +1- Note that cx, c Y , cz, grx, grY, and 
grz were denoted cq^y, cq^x, cxy, 3Q,y, gQ,x, and gxY, respectively, in some previous works. 



APPENDIX C: EXPLICIT CLOCK-FRAME COEFFICIENTS 



This appendix provides the explicit expressions for the clock-frame tilde coefficients in terms of the Sun-frame 
tilde coefficients. For simplicity, we write FIT for the combination f^T and use the abbreviations s x := sinx and 
c x := cos a; for all trigonometric dependences other than the relatively rapid u> s oscillations. 
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The results are as follows: 
b 3 = cosu s T s { ^x(-s Q c c ) + 6y(c Q c c ) + b z (s c ) 



+/3© \d-(-±c a CQC n C£iT 4 
+ b T (-\c a CQC. n cnT + |s q c c sot) + dQ^^CaC^cnr - S(S n c nT - ^SaC^snr) 
+ dxY(-SaC<;c n cnT) + d Y z^c^cnr) + d Z x(-S(s nT ) + gdcaCQC^c^T - SaC^snr) 

C(Sq T ) + H X T(s<;c n cnT ~ c a c c s n cnT) 
+ HYT(-s a c^s. n cnT + s^sot) + HzT(s a C(Cr,cnT - CaC^snr)] j 
| bx(-c a ) + b Y (s a ) 
d-i^SaC^cnr + \c a SQ. T ) + d + (-2s a c n c nT + 2c a s nT ) + b T {\s a c ri cn, T + \c a so, T ) 
+ dQ(ls a c v CQ T - \c a s nT ) + dxYi-CaC^cnr) + d Y z{-s a s v c nT ) + g c (-s a c n c nT - c a s nT ) 
+ gri^SaCrjCQT - feasor) + HxT{s a s v cn T ) + H Y T(~c a s, n c n T) + H Z T(c a c v cnT + SaSnr) j 



sin uoJT< 



+ cos2uj s T s 



|/3 S cL(§c 2a + |c 2 „c 2C ) + 6t(|c 2q + \c 2a C2t) + d Q {j - gc 2 c) + dxv{% s 2a + |s 2q c 2C ) 



+ d Y z(-\c a S2c) + d Z x{\s a S2C,) + .9c(-fc 2ct - \c 2a c 2C ) + g T {- 



\C2o 



d 3 



+ sin2u; s T s {/? s |j_(-±s 2Q c c ) + b T (-^s 2a c ( ) + dxY{\c 2a c c ) + d Y z(^s a s ( ) + <i zx (±c a s c ) 

+ g c {s 2a c c ) + .gT(5 s 2 Q c c )J | 
+ {& d-{-\c 2a + \c 2a c 2C ) + d+{-2) + b T {-\c 2a + |c 2q c 2C ) + d Q {\ - \c 2C ) + d XY {-\s 2a + \s 2a c 2Q ) 

+ dYz(-jc a s 2( ) + d Z x{\s a s 2( ) +g c (\c 2a - ^c 2a c 2C ) + g T (l + \c 2a - ic 2a c 2C )J j, (CI) 

= cosw s T s | d x {~s a CQ) + Jy(c Q c c ) + d z (s^) 

+/3© d-(f C(Sn T ) + d + (-3c a c c c n c nT - 3s ( s v cq T ~ 3s a c ( SQ T ) 

+ b T {-jc a c ( c n c nT - ^s c s n cnT - f s q c c sq T ) + dQ^CaC^cnr + ^s c s n c nT + |s Q c c s n T) 
+ d X Y{\s c^nr) + d-Yzi-s^CrjCnT - ^c a c c s n c nT ) + dzx(s a c<;s v c nT + ^s^sqt) 

+ gT{jC a CQC ri CQ, T + ^s^s^cnr + jS a c^sn T ) + gxY(-^s c c v c nT - CaC^SjjCnr) 
+ gxz(s(C v cn T + ^c a c^s n cn T ) + gYx(-s a c^s n c nT - ^s^sqt) + gYz{\s a cc 3 -n c ^T + s^snr) 
+ gzx(s a c(c n cnT + ^CaC^snr) +gzy(-^ 

+ Hxt^sqc^ccit - \c a CQS v cn T ) + HYT(-ls a c^s v cnT + |s c s f2 T) + H Z T{\s a c^c v CQ,T - IcaC^ssir)] j 
+ sincjsTsj dx(-c a ) + dy(-s a ) 



+ d Q^J SaC-qCnr + jC a s nT ) + d X Y {\c a c r) c^ T - s a s nT ) + d Y zi^SaSrjCnr) + dzx{caSr,cnT) 
+ 9t(—j SaCnCQT + jCaSnr) + gxY(s a s v cnT) + gxz{~\s a s v cnT) + gYx(-c a s v cnT) 
+ 9Yz(^c a s v c nT ) +g Z x(c snr) +gzY{-\ snr) 
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+ H X T{\s a S v CnT) + HyTi-^CaSrjCQT) + HzT{\c a C v CQ.T + 5S Q Sn T )J j 

+ cos2w s T s {/? s d-(-^c 2a - ^c 2a c 2( ) +~b T {jQ - ^c 2C ) + dg(-^ + jqC 2( ) 

+ 4r(-|s2» - jQS 2a c 2<; ) + d~Yz{\c a s 2 Q) + dzx{-\s a S2C,) +9t(~^ + Tg C2 c) 

+ .9xy(|c Q s 2 c) + 9xz{-^c a s 2C ) + gYx(^s a s 2C ) + gYz(-^s a s 2( ) 

+ gzx{-j^s 2a - -^s 2a c 2( ) + gzr(jQS 2a + ^s 2q c 2C )J | 
+ sin2u; s T s {/3 s |j_(§s 2 «c c ) + d X y(-fc 2a c c ) + dyz(-fs a s c ) + <2 zx (-fc a s c ) 

+ gxY{-jS a s c ) + gxz{js a s c ) + gYx{jc a s c ) + gYz{~jc a s c ) 

+ gzx(-jc 2a c c ) + gzY(jc 2a c c ) | 
+ [(3 s [d-(^c 2a - ^c 2a c 2C ) + J+(3) + 6 T (i| - ^c 2C ) + J Q (-i| + j|c 2C ) + d X Y(jQS 2a - ^s 2a c 2C ) 

+ d~Yz(lc a s 2C ) + dzx(-^s a s 2C ) + # T (-±§ + ^c 2C ) + g X Y (|c a s 2 c) + ffxz(-f c Q s 2C ) 

+ <h^(fs a S2c) + .gyz(-|s a s 2C ) +gzx(^s 2Q - ^s 2q c 2C ) + g zy (-^s 2q + ^s 2a c 2C ) |, (C2) 

,g d = cos w s T s | ^£, x (-s a c f ) +g Dy (c a c c ) +5£, z (s c ) 

ccsot) + .g c (2c Q c c c^cnT - 2s a c<;SQT) + gricaCQC^cuT + s^s v cn T ~ SaC^snr) 
+ gxYi-^s^CjjCnr) + .gxz(2c a c c s^co T ) + gYx(-^SQ T ) + gYz{^s a CQS. q cnT) 
+ gzx{2c a CQSQ,T) + ffzy(-2s a CfC I7 cnT)j | 
+ sinw s T s { 5r>x(-c a ) Ijcyf-sJ 

sot) 

+ gxz(~2s a s n cn T ) + gYz(2c a s n cn T ) + gzx(-2s a sn T ) + gzY{-2c a c r) CQ, T ')^ } 
+ COS2WS7; j/? s 6 T (± + \c 2a - \c 2C + \c 2a c 2C ) + .g c (-fc 2a - \c 2a c 2C ) 

+ 9t{-\ - jc 2a + ±c 2C - \c 2a c 2( ) + gxY(^c a s 2( ) + gxz(-\c a s 2C ) + gYx{\s a s 2C ) + gYz{-\s a s 2 c) 

+ gzx{-js 2a - \s 2a c 2( ) +g ZY {jS 2a + js 2a c 2C ) | 
+ sin2w s T s j/3 s [fe T (-s 2 aC C ) +g c (2s 2a c c ) + <? T ( s 2aC C ) + gxY(-s a s ( ) + g X z(saS C ) 

+ ffyx(c Q s c ) + gyz(-c a sc) + gzx(-c 2a c c ) +5zy(c 2a c c )J | 
+ {ft 6t(-| - \c2a - \c 2( + \c 2a c 2( ) + g c {\c 2a - \c 2a c 2C ) +gr(-\ + \c 2a + {c 2C - jC 2a c 2C ) 

+ gxY{\c a s 2C ) +gxz{-\c a s 2C ) +gYx(\s a s 2 c) + gYz{-\s a s 2C ) 

+ gzx{\s 2a - \s 2a c 2C ) +g Z Y{-\s2a + |s 2a c 2C )J |, (C3) 

c q = coscj s r s |/3 s \cTx{2s a cc) + c T y(-2c Q c c ) + c T z(-2s ? )] } 
+ sinw s T 5 {/3 s [cTx(2c a ) + c T y(2s Q ) | 

+ cos2w s T s { [c_(§c 2q + |c 2 „c 2C ) + cq(| - |c 2C ) + c x (-|c Q s 2C ) + cy(|s Q s 2C ) + c z (|s 2a + |s 2a c 2C ) 
+(3® c T x{-^s 2a c n CQ T — ^s 2a c 2 (;C v c nT - jS a s 2 ^s v cnT + | s ot + §c 2q sqt - |c 2 ^sq T + |c 2Q c 2? Sn T ) 



+ c T z(jC a s 2C c v c nT + \s v c nT - jc 2C s v c nT + js a s 2C s nT ) I 



+ sin2w s T s 



+ 



I c_(-§s 2q c c ) + cx(|s Q sc) + cy(§c Q s c ) + c z (§c 2a c c ) 

{ [ S -(-|c2a + l>C 2a C 2C ) + Cq(-| - §C2C) + Cx(-§C a S 2 c) + Cy(|s Q S 2C ) + C Z (-§S 2 a + |s2 Q C 2C ) 
+/3© \c.Tx{^S 2a C v CQ,T ~ ^S 2a C 2Ci C n CQ,T - jS a S2C s v c i1T ~ |so T - ^C 2a S n T - ^C^Sqt + §C2aC2CSOT) 

+ cty(|c^cot - f c 2a c v cnT + ^c 2( ;c v cnT + ^C2 a c 2 (;c v cn T + jc a s 2 QS n Ci-iT — ^s 2a SQ T + ^s 2a c 2 ^sn T ) 

+ C T z{\c a S 2Ci C n CQ,T - \s v C{1T - jC2C s n c ^T + |s q S2C s ot)J j, 



= cos Lu s T s {f3 s [g X Y(s a c c ) + gxz{s a c c ) + gYx(-c a c ( ) + gYz{-c a c ( ) + g zx {-s c ) + <?zy(-s c ) 
+ sinw s T s |/3 s ^xy(c a ) + gxz(c a ) + gyx(s Q ) +gi-z(s Q ) 
+ cos2u; s T s { g_(fc 2Q 



c 2a c 2C ) + .9q(| - §C2c) + .9tx(-|c q s 2< ) +gTy(|s Q S2c) +9Tz{\s 2a + ^s 2a c 2( ) 



+ gc(js a s 2C c v c nT - js 2a s v cnT - f S2aC2CS v ca.T + jc a s2(Sn T ) 

+ gT{\s a S2C,C v CQT ~ §S2aS^CnT ~ ^S 2a C 2 {S n Cot) 

+ CDfY(-f ^aC^C^T - g^a^fC^CfiT — f s f2T + gC2aSfiT + f C2(Sqt + ^C 2a C 2 ^SQ,T) 

+ gxz(-|s Q s 2 ^s^c f2 T + jsqt - |c2csnr) 

+ 9YX{^C V CQ.T + §C2aCr ? CnT - §C2CC^CoT + § C2a C 2 <; C,, Cq T + §S 2q Sqt + ^S 2a C 2 ^Sn T ) 

+ 9Yz(-jC n c nT + jC 2 £C v cnT + jC a S2(S n c nT ) 

+ gzx(^Sr,CnT + ^C 2a S n CnT ~ ^C 2 (S n C nT + \c2 a C2CS v CQ.T + jS a S 2( ;SnT) 

+ gzY(jC a s 2C c, n c nT + ^s n cnT - f c 2a s n c nT - |c 2? s^c t - |c2aC2CS J) Cnr) } 
+ sin2cj s T s | .9-(-§S2 Q c c ) +,g TX (|s Q s c ) +.g T y(§CaSc) +5Tz(|c2aC<;) 

+/3© b T (-%c a S£C n cnT + ^c 2a CQS v Ci-i T ) + g c {^c a SQC ri c n T - 3c 2q c c s, ) csit - ^s a s^sn T ) 

+ gri^CaS^CnT - |c2aC^S r) Cn T ) 

+ gXY{—\ C 2a CQC^C^iT - \s2aCQSQr) + (j X z(-§C a Sf SjjCfir) + ffFX (- §S2aC C C^CoT + ^C2 a C(Sn T ) 

+ 9y zi-^SaS^s^Ciir) + gzx(-^s 2a c i s v cn T + ^c a s c s nT ) + g ZY {-^s a s ( c v CQ T + ^s 2a c c s n c nT ) | 
+ { <7-(-fc2a + lc 2a c 2C ) + ()q{-\ - §C2c) + 9Tx{-jC a s 2C ) + g T Y{jS a s 2C ) + g T z{-ls 2a + |s 2a c 2C ) 
+/3© &t(-§ s a S2(C ri cnT - ^s 2a s, n caT + ^s 2a c 2 QS, q c nT ) 

+ gc(jS a S 2 (;C v CnT + jS2aS v CnT - jS2aC2CS v CnT + jC a S2(SnT) 

+ gT(jS a S2cc n c nT + ^s 2a s v c nT - ^s 2a c 2C s v c nT ) 
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+ gxz{-js a s 2 Qs v c nT - \s nT - jc 2( sn T ) 

+ 9Yx(—^c v cnT — ^c 2a c 71 cnT - §C2cc,,cnT + ^c 2a c 2 QC v cnT — ^s 2a snr + ^s 2oi c 2 (;sqt) 
+ 9yz{\c v cq,t + jc 2 ^c v cq T + jCaS^s^car) 

+ gZx{—\s v CaT - ^C 2a S ri CUT - §C2CSr;CfiT + ^C 2a C2^S v CUT + jS a S 2 £SnT) 



(C5) 
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